Diffractive x Production at the Tevatron and the LHC 
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We present predictions for the diffractive production of \ mesons in the central rapidity region 
usually covered by collider detectors. The predicted cross sections are based on the Bialas-Landshoff 
formalism for both exclusive and inclusive production and makes use of the DPEMC Monte-Carlo 
simulation adapted with kinematics appropriate for small-mass diffractive production. We compare 
generator-level results with a CDF measurement for exclusive \ production, and study background 
and other scenarios including the contribution of inclusive \ production. The results agree with the 
Tevatron data and are extrapolated, highlighting the exclusive \ CQ production at LHC energies. A 
possible new measurement at the Tevatron using the D0 forward detectors is investigated, taking 
advantage of the dominance of exclusive production for high enough diffractive mass fraction. 



I. INTRODUCTION 
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Exclusive and inclusive central diffractiveproduction of heavy states have been studied previously in the double 
Pomeron exchange formalism (DPE) 0, @, i, i, H, H, 0, [1] , and experimental results have been presented Q , attracting 
theoretical attention. Indeed, massive dijets have been copiously produced in inclusive diffractive production at the 
Tevatron and there is presently an active search for exclusive diffractive production of heavy states. One motivation for 
this search is that the Higgs boson could be produced in such a mode, allowing for a good mass determination for this 
' elusive particle if such exclusive events are identified fllj . However the estimated cross sections and signal-over-noise 
ratios are still a matter of debate. 

One way to address this problem is looking for a similar production mechanism with lighter particles like the x 
mesons [10[ . This would give rise to high enough cross sections to check the dynamical mechanisms. Indeed, exclusive 
i-^H production of Xc has been reported by the CDF collaboration [ijj with an upper limit for the cross section of 

•i-H . (Jexc(pp^P+ J/i> + l+p) < 49 ± 18(stat) ± 39(sys) pb . (I) 

X: 

Our goal in the present paper is to analyse both exclusive and inclusive diffractive production of \c and Xb m 
the context of the Tevatron, making useful predictions for the LHC. A possible new measurement using the D0 
forward detectors is also investigated. As a general framework, we use an extension to low mass states of the Bialas- 
Landshoff (BL) model [3 for both exclusive and inclusive production. More precisely, we consider the inclusive model 
as applied in Ref.Q and the extension to purely exclusive processes @, using the corresponding color-singlet (J z — 
0) subprocess cross sections of the original Bialas-Landshoff formalism for diffractive Higgs production [l|, applied 
also to heavy quark pairs in |3j. In fact, in our approach, both inclusive and exclusive diagrams come from the same 
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non-perturbative approach, originated in The model for exclusive production starts with the same soft Pomeron 
exchange diagrams (for ordinary dijet production the g g — > gg diagrams are also included [131 ]) and corrects the 
result by non-perturbative rapidity-gap survival factors [14 , 1 1 51 ] - For inclusive production in the same framework, 
no such factor is applied, but the normalisation is fixed by fitting the CDF dijet data [§]. The energy dependence 
is related to the rise of ordinary hadronic cross-sections through features of the soft Pomeron We think that 
the non-perturbative Bialas-Landshoff approach is well suited for the production of relatively low-mass diffractive x 
states 1 . 

The full kinematics is evaluated for low-mass exclusive production and implemented in the DPEMC Monte Carlo 
generator originally designed for heavy mass states. We use the generator level distributions to estimate the 
measurements at detector level, and to study the systematic uncertainties due the lack of a full detector simulation. 

The paper is divided in the following way: in section [Til we recall briefly the Bialas-Landshoff (BL) formalism for 
both exclusive and inclusive diffractive production and give the theoretical cross section formulae. They include the 
full diffractive kinematics for a low mass state, which is derived in section HTT1 and implemented in the DPEMC Monte 
Carlo. In section [Tvl we show the predictions for the exclusive and inclusive diffractive x production cross sections 
at the Tevatron and LHC. In the subsequent section [V] we compare our predictions with the CDF measurement and 
investigate the possibility of a new measurement at the D0 detector in section IVI1 Section IVIII is devoted to the 
predictions for exclusive Xco production at the LHC and in section IVIIII we present our conclusions and an outlook 
on this interesting mode of production. 



II. THE BIALAS-LANDSHOFF FORMALISM 



One generally considers two types of DPE topologies for the production of a heavy state: exclusive DPE p], 0, H, 0, II] , 
where the central object is produced alone, separated from the outgoing hadrons by rapidity gaps: 



hh — > h + heavy object + h 



(2) 



and inclusive DPE 0, 0, H, d, 0], where the colliding Pomerons are resolved (very much like ordinary hadrons), 
accompanying the central object with Pomeron "remnants" (X,Y): 



hh — > h + X + heavy object + Y + h . 



(3) 



in both cases h represents the colliding hadrons. 

In general, exclusive production is considered most promising, since a better signal-to-background ratio is expected. 
Indeed, if the events are exclusive, i.e., no other particles are produced in addition to the heavy object and the outgoing 
hadrons, the measurement of the scattered hadrons in near-beam detectors provides information on the mass of the 
heavy object [HI], and the dynamics of the hard process. However, the approximations made for heavy mass states 
lead to inconsistent results in the case of lighter objects, so in our case we have to develop the full kinematics for all 
mass states. 

In order to evaluate both inclusive and exclusive diffractive production, we use an extension [8[ to the purely 
exclusive processes of the original Bialas-Landshoff formalism for diffractive Higgs production [l|, also applied to 
heavy quark pairs Q . In this extension, both inclusive and exclusive diagrams come from the same approach and are 
based on soft Pomeron exchange diagrams (for ordinary dijet production the gg — > gg diagrams are also included [Hi]). 
In the exclusive case one has to correct the result by non-perturbative rapidity-gap survival factors [3, EH , while 
in the inclusive case the normalisation is fixed [H] by reference with dijetproduction @. The values of the survival 
probabilities are the same as for the diffractive Higgs production in Ref. [8|, namely 0.1 and 0.03 for respectively the 
Tevatron and the LHC. The energy dependence of the cross section is related to the rise of ordinary hadronic cross 
sections through features of the soft Pomeron [l6| • 

With respect to previous works, the formulae of the Bialas-Landshoff cross sections have to include now the full 
kinematics for diffractive x meson states and thus valid for small masses in general. This kinematics is explained in 
the next section. For the exclusive x meson production in hadron-hadron collisions, one has 
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(4) 



At this point, only XcO and XbO mesons can be evaluated. Indeed, due to non-zero p± of the initial state (pp or pp), there will be 
contributions of \c 2 an d Xb 2 m the exclusive amplitude. This will be discussed in a further paper. 



where C x is the normalisation constant (cf. [l| for the Higgs boson) , My is the \ meson mass and A x the slope at the 
Pomeron vertex, a 1 and e are the standard soft Pomeron parameters [161 ] . The gluon-gluon \c coupling satisfying the 
J z — rule is taken as in Ref. fiol ] . The parameter C x contains a non perturbative part due to the gluon coupling G 
to the proton. We kept the original value in the model of Bialas Landshoff which means G 2 /An ~ 1. This is of course 
an order of magnitude, and it leads to an uncertainty on the exclusive production cross section. 

Most importantly, M^iff is the expression of the diffractive mass produced in the central region in terms of the 
kinematic variables. For low-mass states, as the x's, it sensibly differs from the expression for heavy diffractive states, 
as will be derived in the next section and in Appendix [A"l This is one major modification one has to introduce w.r.t. 
the original formalism. 

For inclusive production, the cross sections are given by 

/ T S T 9\2e / M 2 \ 



where one makes use of the parton structure functions G p {x 9 i ,iJ,) in the Pomeron, for a given scale /j,. Note that in 
Equation (JSJ), the normalisation C™ c is fixed by normalising the prediction on the measurement of the dijet diffractive 
cross section by the CDF collaboration 

In both equations, the variables Vi and denote the transverse momenta and fractional momentum losses of the 
outgoing hadrons. In the second equation, xf denote the fractional momentum carried by the gluons in the Pomeron. 

Note that in the inclusive case, the usual kinematics is used for e 



III. FULL KINEMATICS FOR EXCLUSIVE PRODUCTION 

Let us first recall the method used to generate exclusive events in the DPEMC Monte Carlo generator [ItJ • The first 
step is to randomly generate t\, t% and £i, following an exponential distribution for |tj|, where |tj| is the 4-momentum 
transferred and £j is the momentum loss for the hadron i. Exclusive events have the property that the full energy 
available in the center-of-mass is used to produce the diffractive object, or in other words there is no Pomeron remnant. 
The diffractive mass can be expressed as: 



(6) 

The value of £2 is thus imposed by this relationship. The produced events in the generator are then weighted as 
usual by the cross section. 

The approximation ^ is no longer true for low mass states such as \ mesons, and we had to modify the method to 
generate events in this case. Wc derived the diffractive mass from full 4-momentum conservation (see Appendix A). 
Using the full kinematics, Equation © is replaced by 



m%„ _ . x (1 + t' -«f -<■> (, _ n, _ + !_&)) , (7) 

11 \ ZCOSU1COSO2 \ COSt>l COS02 J J 

where 57 = —cosQ\cosQ2 + sind\sinQi{cosLp\cos<p2 + simpisin^), is the scattering angle and (p the polar angle. It 
is important to notice that this formula depends not only on £1 and £2 but also on the angles of the hadrons #i,</>i 
and 02,<f2- t and are related by the following formula: 



We use the following method to generate low mass exclusive events within the framework of DPMEC [TtJ • We start 
by generating 9i, 62 (following an exponential distribution) and £1 randomly, which gives t\ by the equation ([S]). £2 
is then computed inverting equation ([7]), giving also £2- The events are then weighted according to the cross section. 
The new steps are thus to use the variables 8 and £, avoiding the cumbersome solution of equation ([7|) in terms of £ 
and t. 



IV. EXCLUSIVE AND INCLUSIVE X c , Xb PRODUCTION CROSS SECTIONS 



Table (p} presents our results for the cross section predictions at the Tevatron and the LHC. The energy dependence 
for the exclusive production cross section can be seen in the Fig. [T]for the Xb (left) and Xc (right) mesons. The gap 
survival probability (the probability of the gaps not to be populated) Sg ap is taken to be 0.1 at the Tevatron and 0.03 
at the LHC (cf. A.B. Kaidalov et al in [3). 



<j(nb) 


Tevatron y/s = 1.96 TeV 


LHC y/s = 14 TeV 


cr exc(Xco ) 


1.17xl0 3 


0.804 xlO 3 


<Texc{Xb ) 


4.4 


3.29 


cr inc(Xco) 


1.8 xlO 4 


4.8 xlO 4 


ow(xf>o) 


20 


1.8 xlO 2 



TABLE I: Cross sections (in rib) for exclusive and inclusive production at the Tevatron and the LHC. 



Note that our estimates for exclusive X c production at the Tevatron are higher than other recent predictions [1 01 ] - 
The differences arise mainly because of the distinct Pomeron fluxes in the models, which means that the BL exclusive 
cross section has a smoother dependence with the center of mass energy (\/s) (see Fig. [T]). We should also note that 
our model does not contain Sudakov factors contrary to Ref. [J, [1(3] . The effect of the Sudakov suppression is however 
supposed to be small at small masses, and plays a more important role at higher masses. 




FIG. 1: Exclusive Production cross section of Xbo (left) and XcO mesons (right). 

Possible sources of background for exclusive production include: 

• cosmic events which can fake an exclusiveness; 

• exclusive and inclusive X i production; 

• quasi-exclusive x production, which are defined by a high mass fraction Fm (e.g. Fm > 0.75-0.95), 

where the mass fraction Fm is the ratio between M x and the total diffractive mass Mdiff. We will study this last 
contamination to the exclusive signal in the following. 



V. COMPARISON WITH THE CDF LIMITS ON X c PRODUCTION AT THE TEVATRON 



As shown in table HI exclusive Xc production at the Tevatron has a noticeable cross section, with the potential to 
be detected in the traditional decay Xc —* J/i/)(—> H + + fJ>~) + 7- The CDF Collaboration has presented preliminary 
results [H[ for exclusive J/ip + l production using the rapidity gap selection of diffractive events in Run II (y/s =1.96 
TeV). The cuts used by CDF are the following: pt(a* ± ) > 1-5 (GeV), [77(7) | < 3.5 and |r](^ ± )| < 0.6. The CDF 
collaboration provided an upper limit on \c a production cross section at the Tevatron, uncorrected for residual 
backgrounds such as cosmics 2 : 

VexcipP -> P + J/tp + 7 + P) < 49 ± lS(stat) ± 39(sys) ph. (9) 
If we apply the CDF cuts at generator level, we predict the following cross section 

Oexc{PP -> P + Xca{~* J/ih) +P) = 61p6. (10) 

In order to make a more realistic comparison, we need also to consider the non-exclusive background which can 
enter directly in the experimental cross section determination. In particular, the contamination due to quasi-exclusive 
events need to be considered properly. In this class of events the QCD radiation, or in other words the energy loss 
in the Pomeron remnant and in soft QCD radiation, is small. CDF remove most of the inclusive background using a 
cut on the mass fraction, Fm > 0.85. For exclusive events, Fm is expected to be close to one, since the full available 
energy is used to produce the \ meson, and smaller than one for inclusive events. However, this cut is of course 
applied at the detector level by the CDF collaboration, whereas we can only apply it at the generator level. Due to 
the fact that we are missing the smearing between detector and generator levels, we choose to investigate the effect 
on the cross section due to various mass-fraction cuts, as displayed in Table [TT1 

We also compared the transverse momentum (pr) distributions of the \ Co for inclusive, exclusive and quasi-exclusive 
production, the latest defined as Fm > 0.75 (Fig. [2]). As expected, the exclusive pt distribution reaches higher values, 
so this variable could be used to enhance the exclusive production signal. 




FIG. 2: Transverse momentum of XcO for exclusive production (solid line), inclusive production (dashed line) and for quasi- 
inclusive events (filled) at yfs = 1.96 TeV (the normalization used a luminosity of 100 nb^ 1 ). 

The quasi-exclusive contamination to the exclusive signal also suffers from the unknown gluon density in the 
Pomeron at high f3 in particular. This uncertainty, estimated to be about 50% (l8j . can be taken into account by 
multiplying the gluon density in the Pomeron, measured at HERA, by a factor (1 — [3) v where v varies between -1.0 
and 1.0 [19(. If v is negative, we enhance the gluon density at high (3 by definition. The QCD fits to the HERA data 
lead to the following value of the v parameter: v — 0.0 ± 0.6. 



2 This cross section has been corrected for the J/^ branching ratio into muons. 



In summary, the sources of contamination to the exclusive signal can derive from the smearing in the mass fraction, 
which is a pure experimental effect, or the uncertainties in the gluon density in the Pomeron, which is a theoretical 
limitation. 

Table II gives the quasi-exclusive cross section at the Tevatron after applying the CDF cuts, and different cuts on 
the mass fraction and values of the v parameer. We thus find that the signal seen by the CDF collaboration could be 
explained by a combination of a higher gluon density at high f3 and some smearing effects due to the reconstruction 
of the mass fraction. Table III gives the exclusive production cross section after each of the CDF cuts. 



Mass Fraction Cut 


v = 


v = -1 


v = -0.5 


v = 0.5 


v = 1 


> 0.75 


14.33 


194.94 


52.28 


3.88 


0.84 


> 0.8 


5.4 


118.87 


27.15 


0.84 


0.17 


> 0.85 


2.02 


61.89 


11.13 


0.17 





> 0.9 


0.34 


28.43 


2.87 








> 0.95 


0.08 


19.48 


0.84 









TABLE II: Quasi-exclusive cross section (in pb) at the Tevatron, after CDF cuts, using different Fm and gluon distributions. 



CDF cut 


1 


2 


3 


4 


5 


Exclusive cross section (pb) 


5.56 xlO 3 


7.97 xlO 2 


5.25 xlO 2 


61.47 


61.21 



TABLE III: Exclusive cross section a exc (pp — * p + \co{~* J/^l) + P) (in pb) at the Tevatron energies for each CDF cut: 1 - 
one muon with pr > 1.5; 2 - one muon with pr > 1.5 and \r/\ < 0.6; 3 - two muons with pr > 1.5; 4 - two muons with pr > 1.5 
and \rj\ < 0.6. 5 - same constraint of the forth column plus one gamma with 77 < 3.5. 



VI. POSSIBILITY OF A NEW MEASUREMENT AT D0 

We now examine the possibility of measuring the exclusive Xc production at the Tevatron using the roman pot 
detectors in the D0 collaboration. 

After interfacing the DPEMC generator with a program designed to propagate (anti)protons through the Tevatron 
lattice, from the D0 interacting point to the roman pot detectors [20I ] , we can predict the number of events observed 
for different tagging configurations. The Forward Proton Detector (FPD) installed by the D0 collaboration [2l| 
consists of eight quadrupole spectrometers, four being located on the outgoing proton side, and the other four on the 
antiproton side. 

Each spectrometer allows the trajectory reconstruction of the outgoing protons and antiprotons near the beam pipe 
determining their energies and scattering angles. The quadrupole detectors are sensitive to outgoing particles with 
\t\ > 0.6 GeV 2 and £ < 3.10 -2 , with good acceptance for high mass objects produced diffractively in the D0 main 
detector. 

We use the following selection cuts: ((pt(m + ) > 2.0 (GeV) or p T (n~) > 2.0 (GeV)) and |t7(^ ± )| < 2.0 and 
|»7(7)| < 3.0) (see Table US}. 



Regular Tevatron Stores - L = lOOpb 1 


Scenario 


A 


B 


C 


D 




D0 selection 


1.2 xlO 8 
1.8 xlO 2 


2.6 xlO 6 

2.7 xlO 1 


4.8 xlO 6 
3.0 xlO 1 


2.9 xlO 5 
1.5 



TABLE IV: Number of exclusive Xc events at the Tevatron (MC error ~ 10% ) for a regular Tevatron store. The scenario 
represents all decay channels included without selection cuts. The columns represents the number of events: A - all (without 
p or p tagging); B - tagged in the p side quadrupole; C - tagged in the p side quadrupole and D - double tagged events in the 
quadrupoles. 

We note that the number of events in double tagged configuration is quite small after applying the selection cuts, 
so this configuration might not be useful. However, a single tag event with a rapidity gap on the other side yields 
a good number of events, with the additional benefit of having the kinematics determined for one of the scattered 
particles. 



VII. EXCLUSIVE xc PRODUCTION AT THE LHC 



From the results on section IV, we expect a high production cross section of \c mesons at the LHC. We can estimate 
the number of events accessible to the TOTEM/CMS detectors, benefitting from the good acceptance of the TOTEM 
detectors for low mass objects. The TOTEM acceptance for the high (3* optics and low £ values is typically 90 %, 
for the range < \t\ < 1 GeV 2 . Then for 10 pb -1 of data, 5.3xl0 6 double tagged events are predicted, with no 
requirement in the central detector activity. In this way, one might look for the \c a in the reconstructed diffractive 
mass. 

If central activity is required, the lowest possible muon px cut at low luminosity is on the order of pt > 1.5 
(GeV) for |r?| < 2.4 Otherwise, the muon pt threshold would be 4 (GeV) , which is too high for exclusive \ CQ 

production. The predictions for exclusive and quasi-exclusive production at the LHC are shown in tables [V] and I VII 
We note that the number of events can be dominated by exclusive production, independent of uncertainties in the 
gluon distribution, if a high enough cut on the mass fraction can be made, for instance at 0.95, which requires a good 
coverage of the CMS detector at high rapidities. 



Mass Fraction Cut 


v = 


v = -1 


v = -0.5 


v = 0.5 


v = 1 


> 0.9 

> 0.95 


1.35 



138.11 
13.83 


17.88 
1.18 


0.34 



0.17 




TABLE V: Quasi-exclusive cross section (in pb) at the LHC, after central activity cuts, using different mass fractions and gluon 
distributions, defined in section V. 



Central cut 


1 


2 


3 


4 


Total 

After Totem Acceptance 


3.74 xlO 3 
3.03 xlO 3 


1.43 xlO 3 
1.16 xlO 3 


3.64 xlO 2 
2.95 xlO 2 


1.27 xlO 2 
1.03 xlO 2 



TABLE VI: Exclusive cross section (in pb) a exc (pp — » p + \ ca {— » J/ip^f) + p) at the LHC energies for each central cut: 1 - one 
muon with pr > 1.5; 2 - one muon with pr > 1.5 and \r]\ < 2.4; 3 - two muons with pr > 1.5; 4 - two muons with px > 1.5 
and < 2.4. 



VIII. CONCLUSION 

We calculate the diffractive production cross section for \ mesons at the Tevatron and LHC using an extended 
version of the Bialas-Landshoff model, including the full kinematics needed for low mass states. Both exclusive and 
inclusive production have been evaluated and discussed with various physically motivated cuts at generator level. The 
results can be listed as follows. 

• The results for exclusive production at the Tevatron agree with a recent CDF upper limit for the exclusive 
production of Xc , with the default parameters of the model (assuming that the non perturbative gluon coupling 
to the proton G 2 /47r ~ 1); 

• In the same conditions, the non-exclusive background (in particular "quasi-exclusive" events which we evaluate 
from the inclusive spectrum) can reach similar levels as the exclusive signal, due to experimental and theoretical 
uncertainties, making difficult to establish the observation of exclusive events within the current parameters; 

• The possibility of observing exclusive Xc production at the Tevatron, using the D0 forward detector is thor- 
oughly investigated, showing the possibility of a measurement if a tight cut on the ratio between the \c mass 
and the total diffractive mass can be performed successfully. 

• Exclusive production at the LHC, using the CMS/TOTEM detectors, is also investigated and appears promising, 
but again a high enough cut on the mass fraction is to be made before the observation of exclusive production 
can be established with a good confidence level. It could be suitable for a mass fraction, for instance at 0.95, 
which requires a good coverage of the CMS detector at high rapidities. 
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APPENDIX A: EXCLUSIVE PRODUCTION KINEMATICS 



Due the low mass of the x mesons, we must evaluate the full kinematics of an exclusive production. The start point 
is the 4-momentum conservation in the center-of-momentum frame. Using the approximation m p = nip = 0, where p 
and p represent the colliding particles, we can make E p p = \kp p\, and show that: 

s = (k p + kpf - {k p + kp) 2 + M% iff + 2E M (k p + k p ) - 2k M ■ (k p + kp) (Al) 

where k Pt p is the particle 3-momentum, s is the center-of-mass energy, Mdiff is the diffractive mass, and Em and ku 
are the respective energy and momentum of the produced state. 

Due to energy and 3-momentum conservation Em = \fs ~ k p — kp and Um — ~(k p + kp). Moreover, we define 

kp ' kp = kp kp ^ A. 2 ^ 

where = —cos9 p cos9p + sin9pSin9p(cosip p cosipp + sinip p siri(p p ), 9 is the scattering angle and (p the polar angle. 
Using equation (|A2[) and the conservation constraints in equation (|A1[) . it can be shown that 



Thus 



using the definition of £ 



2k p kp(l ~n) + Mj iff + 2^(kp + k p ) - Akpkp{l - fi) (A3) 
M hf = s + 2k p kp{l -CI)- 2«/s(h p + k p ) (A4) 



£p,p - 1 kl nitial ^ kp,p ~ cos9 p p^ 1 ^ A5 ^ 



Using equation (|A5|) in (|A4| : 



s 2cos9pCOSp \ cos9 p cos9p J 



In the case of |i Pl p| — * 0, which means 9 Py p — > 0, the relation M^** — s£ p £p is obtained. 

Note that for convenience in the main text, the proton kinematics are labeled by the index 1 and the antiproton by 2. 
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